Adiabatic quantum computation with Rydberg-dressed atoms 
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We study an architecture for implementing adiabatic quantum computation with trapped neutral atoms. 
Ground state atoms are dressed by laser fields in a manner conditional on the Rydberg blockade mecha- 
nism, thereby providing the requisite entangling interactions. As a benchmark we study the performance of 
a Quadratic Unconstrained Binary Optimization (QUBO) problem whose solution is found in the ground state 
spin configuration of an Ising-like model. We model a realistic architecture, including details of the atomic 
implementation, with qubits encoded into the clock states of 133 Cs, effective B-fields implemented through 
stimulated Raman transitions, and atom-atom coupling achieved by excitation to the IOOP3/2 Rydberg level. 
Including the fundamental effects of photon scattering, we find the fidelity of two-qubit implementation to be 
on the order of 0.99, with higher fidelities possible with improved laser sources. 



In this letter, we consider a new platform for adiabatic quan- 
tum computation (AQC) [1| based on trapped neutral atoms 
whose coupling is mediated by the dipole-dipole interactions 
of Rydberg states. The algorithm is implemented by contin- 
uous transformation of the Hamiltonian from an initial form 
whose ground state is easy to prepare to the final form whose 
ground state encodes the output of the algorithm. If the energy 
gap between the ground and excited states is sufficiently large, 
the transition from initial to final Hamiltonian can be accom- 
plished efficiently. AQC is particularly attractive because the 
existence of an energy gap can make the system inherently 
robust to certain types of errors. 

In contrast to quantum circuit implementations where 
atoms are excited to the Rydberg state with a resonant 7t- 
pulse [2 -4 1, here we base our proposal on off -resonant Ry- 
dberg dressing of the atomic ground state, studied previously 
in the context of dipolar gases 016]. This leads to an entan- 
gling mechanism that is more compatible with AQC, where 
interactions are always on and can be continuously changed to 
transfer from the initial to final Hamiltonian. Such adiabatic 
evolution has been employed in recent cold atom/ion exper- 
iments to study quantum simulations of Ising models Q [8] . 
As a specific example, we will show how our architecture can 
be used to implement the quadratic unconstrained binary op- 
timization (QUBO) algorithm. We will model the physics of 
its implementation to benchmark the performance of a proof- 
of -principle realization for a few qubits with nearest-neighbor 
interactions. 

The goal of the QUBO problem is to find the A^-tuple of bi- 
nary variables, x = (x\,Xz, ■ ■ • ,xn), Xi 6 {0, 1}, that minimizes 
the function 
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This is equivalent to solving for the ground state of a 
generic Ising model, a problem that is generally NP-hard J5J- 
Nonetheless, specific instantiations of this problem map onto 
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a variety of satisfiability and related algorithms which are 
tractable, and thus provide useful testbeds for the AQC ar- 
chitecture ifTUl . To map QUBO onto the quantum problem, 
each binary variable is replaced by a projector acting on a 
qubit, Xi =>■ (I + a®)/2, where the Pauli matrices are defined 
as usual on the qubit pseudospin, |0) = | t)> |1) = | V)- The 
solution to QUBO maps onto finding the ground state of the 
"problem Hamiltonian," Hp, in the Ising form 

H P = f t haf + £ JijG® ® ai j) , (2) 
''=1 <,;'=i 

where Ty = Jij/4- and hi = hj/2 + Since the Hamil- 

tonian commutes with all ol'\ the ground state is one of the 
computational basis states, which can be read out directly. 

To implement AQC, we reach the ground state of Hp 
through adiabatic evolution. The beginning Hamiltonian, Hp, 
must have a simply preparable ground state with [Hp,H p ] 7^ 0. 

We choose Hp = — , whose ground state is all spins 

up along the ^-direction of the Bloch sphere, |+)®", with 
|+) = (|0) + |1)) /\2. In a linear ramp, we transition from 
Hp to Hp, H{t) = (l-t/T)H B + {t/T)H P , with the time T 
chosen sufficiently long to ensure adiabatic evolution. Opti- 
mized, non-linear ramps between Hp and Hp can improve the 
fidelity of the adiabatic evolution, but the nature of such ramps 
will depend on the problem at hand. 

As a benchmark for performance of this architecture, we 
will study a class of QUBO problems corresponding to 
a one-dimensional spin chain with symmetric interactions, 
J<ij> = J, where <ij> denotes nearest neighbors. We choose 
the values /z; to be equally spaced and less than J, hi = i?>E 
with N8E < J for N qubits. The solution to this QUBO 
problem is trivial; minimization is achieved with the state 
1 1010 • • • 10) for even N or |0101 •10) for odd N, i.e., the 
bits alternate between 1 and and the final bit is 0. Fur- 
ther, the gap between ground and first excited states scales as 
N~ l , so the necessary evolution time to maintain adiabaticity 
grows linearly. We consider this example only as a proof-of- 
principle of the method that can be modeled numerically for a 
few qubits. 
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To implement this test-bed algorithm in a neutral-atom sys- 
tem, we consider cesium atoms that are optically trapped 
along a line. Qubits can be encoded in two hyperfine magnetic 
sublevels in the ground-electronic state of alkali-metal atoms, 
e. g., the "clock states" of 6S 1/2 133 Cs: |0) = \F = 4,M F =0), 
|1) = \F = 3,Mf = 0). Arbitrary single qubit Hamiltonians 
of the form H = B • a can be achieved with stimulated two- 
photon Raman transitions in the standard manner, with negli- 
gible photon scattering over the duration of the evolution for 
sufficient detuning and intensity of the lasers. The last critical 
ingredient is the coupling matrix of pairwise interactions, Ju. 

To generate interactions between widely separated neutral 
atoms, we dress the ground state using the Rydberg blockade 
mechanism as discussed by Johnson and Rolston 0. The key 
idea is to induce a light shift (LS) on the atoms that depends 
on the dipole blockade. Consider a simple model of two-level 
atoms, with ground-state \g) and excited Rydberg state \r), 
and for interatomic separations such that the Rydberg block- 
ade is perfect. Double occupation of \rr) is suppressed, and 
the laser couples the ground \gg) state to the excited "bright 
state" (\rg) +e ,A *|gr)J /\/2, an entangled state, where A(j) is 
the phase difference of the lasers at the positions of the atoms. 
The residual difference between the LS of blockaded and non- 
interacting atomic pairs is the desired coupling constant, 
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Here 2E^ is the noninteracting component of the two-atom 
LS, which can be removed via single-atom addressed Raman 
lasers. Note, for a perfect blockade, the interaction strength is 
independent of the motional phase, e' At>f , as long as the detun- 
ing is large compared to the Doppler width. 

For application to quantum computing, we require this in- 
teraction to be conditional on the states of the qubits. For 
the Ising (QUBO) problem, we seek to implement the pair- 
wise coupling fijoi'^ <X> . For two neighboring atoms, 
we can achieve this when the detuning of the Rydberg laser 
is small compared with the ground-state hyperfine splitting 
(9.2 GHz for 133 Cs) so that the LS is negligible for all but 
a given computational state of the targeted qubits, |xi*2), 
with x E {0,1}. As above, this two-atom ground state 
is dressed through off-resonant coupling to the bright state 
{\ rx 2} +e' A ^\xir)) /a/2, and the effective interaction Hamil- 
tonian is H\ nt w J\x\xi) (x\xt\. Up to single qubit terms (that 
can be compensated by individually addressed atomic LS), 
Hi nl ±(J/4)c z ®G z . The positive/antiferromagnetic (nega- 
tive/ferromagnetic) sign is achieved when x\ — X2 {x\ 7^X2). 
The ability to choose the signs of the elements of Jy provides 
extra flexibility in this platform, even if the sign of the phys- 
ical coupling is fixed in the dressing interaction. These con- 
cepts are illustrated in Fig. 1 . 

One fundamental limitation on the fidelity of operation is 
the scattering of photons due to excitation of the Rydberg 
state at a rate y r = N r T r , where N r is the population in the 
Rydberg state and F r /2n is its linewidth. Scattering of trap- 
laser photons can be made negligible with blue detuned "bot- 
tle traps" iflTI . While F r °c «~ 3 points to larger principal quan- 



tum numbers, a variety of practical considerations limits the 
value of n, including the linewidth of the Rydberg excitation 
laser, the sensitivity of the Rydberg state to ambient fields, 
and the sheer size of the Rydberg atom. For example, for the 
100P state, the quadratic Stark shift due to the weak DC elec- 
tric field is approximately —98 GHz/(V/cm) 2 , and the radius 
of the atom r « 0.7 /jm. This sets a reasonable operating point 
for our architecture. To illustrate our protocol, we consider 
here the 100P 3/2 state, for which T r /2% = 530 Hz. By di- 
rectly dressing the ground state with the Rydberg state using a 
single optical field at A, = 318 nm, we avoid the strong photon 
scattering that arises in the conventional two-photon excita- 
tion scheme via an intermediate excited state and reduce the 
total photon scattering rate by a factor of 10 or more. 

A second fundamental limitation is the accuracy with which 
we can implement the QUBO Hamiltonian using the Rydberg- 
blockade interaction. Ideally, we would like to introduce only 
the desired pairwise couplings specified by the matrix 7y. In 
practice there will be additional perturbations due to the long- 
range nature of the dipole-dipole interaction and the strong 
blockade mechanism. For our geometry, this means that there 
are residual next-nearest neighbor couplings and many-body 
effects (e.g., 0~2|) when more than two atoms are close to the 
blockade radius. Both interaction types will add unwanted 
terms to our final Hamiltonian, potentially shrinking the min- 
imum gap or even changing the final ground state if they are 
too large. However, as long as these effects can be treated as 
a perturbation that is sufficiently small compared to the min- 
imum energy gap, they will not interfere with the adiabatic- 
ity of evolution, and the algorithm will still give the correct 
answer; this sets a minimum acceptable energy gap and, by 
extension, constrains the size of problem that can be solved. 

We must address further details in order to understand 
the exact nature of the coupling J. The qualitative discus- 
sion above holds only for a simplified model of participat- 
ing atomic levels and for a perfect dipole-blockade, i.e., when 
the probability of simultaneously exciting two adjacent Ry- 
dberg atoms is zero. To obtain a more accurate description 
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Figure 1. (Color Online) (a) Schematic of relevant energy levels and 
laser couplings in cesium atoms: Qubits are encoded in the hyperfine 
clock states, controlled by two-photon Raman lasers. Interactions 
between qubits are mediated by off-resonant excitation near a highly 
excited Rydberg state |100P//), tuned to dress either 1 0) or 1 1 ) . (b) 
Dressed Rydberg interaction for two atoms: The logical state |jtiJt2), 
16 {0, 1}, is dressed by the bright state, with one atom in the Ry- 
dberg level; the doubly excited state is blockaded. The result is a 
two-atom light shift with an entangling component J. 
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Figure 2. (a) Doubly excited levels as a function of distance between 
atoms that asymptote near the atomic pair IIOOP3/21 lOOi^rt}- (b) 
Weighting of the levels in (a) by their oscillator strengths to dress the 
two-atom ground state. 



we can find the dressed-state eigenvalues by diagonalizing the 
two-atom system in the presence of the laser field, yielding 
a position dependent J(r) 0. Outside the blockade radius, 
the result is J(r) °< r~ k , where k = 3 for the Forster regime 
or k = 6 for the Van der Waals regime. As we are consider- 
ing direct excitation to a p-state, there may be concern that 
pairs of atoms would couple to noninteracting "Forster zero 
states" that evade the Rydberg blockade 1131 , Such zeros are 
avoided, however, in a more complete description of the elec- 
tric dipole-dipole interaction (EDDI) since mixing occurs not 
only between p- and s-states but also with nearby ti-states and 
higher angular momentum orbitals. 

Inside the so-called blockade radius the situation becomes 
significantly more complex. Our description has been based 
on weak EDDI mixing of nearby Rydberg levels. Under this 
model, inside the blockade radius, the interaction strength 
J(r) decreases very slowly with r, plateauing to the value 
given in Eq. <[3j when the blockade is considered perfect J5). 
For short separation distances, however, this model breaks 
down, as the magnitude of the EDDI grows large relative to 
the bare atomic level splitting. The result is an admixture of 
many more nearby Rydberg states with different nl quantum 
numbers, leading to a splitting of the excited states into a large 
"spaghetti" of energy levels that have more "molecular" than 
"atomic" character fl4*l . 

Figure|2^ shows the portion of the doubly-excited spectrum 
that asymptotes at large separation near to the energy of two- 
atoms in the IOQP3/2 state. This spectrum is calculated using 
a basis of two-atom Rydberg levels |«aWa'm/a) \ n B^BjB m iB), 
with 96 < n < 101, < Z < 6, and including fine structure 
and all magnetic sublevels, a total of ~ 45,000 bare states. 
Clearly, without proper care, this dense fan of excited levels 
can shift into resonance with the laser field, potentially negat- 
ing the blockade effect at the heart of our protocol. The impact 
of such resonances on the dressed states will depend on the os- 
cillator strengths for the transition that connects the two-atom 
ground state to these doubly-exited states via the intermedi- 
ate singly-excited bright states. Figure |2j) shows a plot of the 
spectrum, weighting each level by this oscillator strength. We 
see that outside r«8 ym there is little oscillator strength for 
blue detunings on the scale of interest. Inside this radius, the 
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Figure 3. (Color online) Ground-state interaction strength J(r) as a 
function of distance between the two cesium atoms in F = 4 clock 
states, for Cl r = 10 MHz and A,. = 8 MHz. Calculated using all I < 6 
atomic orbitals (s,p,d,f,g,h,i; purple solid line) and more approxi- 
mately using only the nearest I < 2 orbitals (s,p,d; blue dashed line). 
As seen in the more exact calculation, below r = 8 /jm there are res- 
onances whose exact positions cannot be predicted without taking 
higher-/ states into account. 



situation is much less clear. While there may be regions of 
small oscillator strength for shorter interatomic distances, our 
calculations based on the bare atomic basis fails to converge, 
and the results cannot be trusted. We will restrict our attention 
here to r > 8 fjm, which will give us sufficient coupling and 
control over the atoms. 

Given the doubly-excited spectrum, we can calculate the 
dressed-ground-state coupling J(r). An important considera- 
tion is the choice of detuning of the 318 nm Rydberg laser. 
We seek to maximize the figure of merit K — J/y r so that we 
achieve a large gap between the ground state and excited com- 
putational states of the problem Hamiltonian, while minimiz- 
ing photon scattering over the duration of the evolution. In the 
simplified model that leads to Eq. j3J, |/| continues to grow 
as |A| decreases, while y, saturates at half of the spontaneous 
emission rate. This suggests that the best operating point is 
not one with weak dressing, but with a strong admixture of 
Rydberg character in the dressed-ground state that arises for 
small detuning. For such strong dressing, the minimum possi- 
ble detuning is set by the requirement that the dressed-excited- 
Rydberg state is sufficiently gapped from the dressed-ground- 
electronic states, so that adiabatic evolution can proceed. 

Including the full doubly-excited spectrum shown in Fig. [2] 
the choice of detuning is found empirically. For a Rydberg 
laser that achieves a Rabi frequency Q. r = 10 MHz, we find 
that a good choice of detuning is A r = 8 MHz. Figure [3] 
shows a calculation of J(r) for these parameters, and its com- 
parison to the simplified two-level atomic model. For tightly 
trapped separated atoms, J(r = 8 fjm)/2K — —470 kHz. At 
such a laser power and detuning, there is substantial dress- 
ing, with as much as ~ 20% of Rydberg character in the 
dressed ground states. The maximum photon scattering rate 
is Yr/27t ~ 100 Hz, yielding an excellent figure of merit for 
AQC. Next-nearest-neighbor and three-body interactions for 
these parameters are smaller than the minimum gap for up to 
five atoms; increasing A r and r could increase the maximum 
problem size farther at the expense of K. 
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With these parameters we model the performance of proof- 
of-principle experiments to solve the QUBO problems via 
AQC. The basic protocol is as follows: We optically pump the 
atoms into a clock state and initialize the qubits in an eigen- 
state of o A through the application of a Raman-resonant 7t/2- 
pulse. We then phase shift the Raman beam by Jl/2, leaving 
the atoms in the ground state of the beginning Hamiltonian, 

Hb = —B x Y.i<5x\ where 2B X is the Raman-Rabi frequency. 
This initializes the quantum register in an equal superposition 
of all computational basis states. The transition from initial 
to final QUBO Hamiltonian is achieved by ramping down the 
Raman laser power while ramping up the individual atoms' 
Raman detunings that create the local Hamiltonians hi<J^\ Si- 
multaneously, we linearly increase the Rydberg laser power 

that creates the coupling Hamiltonians JjjG^ <8> o% /4. At the 
end of the ramp, the answer to the algorithm can be read out 
using state-dependent resonance fluorescence. 

The coupling matrix is induced by the dipole blockade 
mechanism discussed above. For QUBO problems under con- 
sideration, we take /y = J V;=;±i> achieved when all atoms 
are arranged in an evenly spaced line, with spacing between 
nearest neighbors d such that J(r) ^> J(2r). As can be seen 
in Fig. [3] r«8 /jm for blue detuning fits this criterion. Note, 
since the coupling parameter /<;/> is positive, while the phys- 
ical J is negative for our system, we achieve the desired an- 
tiferromagnetic Ising coupling by using Rydberg laser fields 
that individually address the atoms, alternately dressing near- 
est neighbors in |0) and 

We take as our parameters B x = /<//> = 470 kHz, and 
hi = (i/N) 118.5 kHz for N qubits, achievable with the atom- 
laser interactions discussed above. The ramp time is taken to 
be 35 /js, sufficiently long to maintain adiabatic evolution, but 
sufficiently short compared to the photon scattering time. We 
treat spontaneous emission from the Rydberg level as effec- 
tively randomizing the magnetic spin state as the population 
cascades back to the ground electronic ground state. For prac- 
tical reasons, the detection scheme does not distinguish be- 
tween different magnetic sublevels in the same hyperfine sub- 
space. All magnetic sublevels in F = 3 are treated as logical- 



1 and those in F — 4 as logical-0. Our simulation for two 
qubits, with the correct solution to QUBO encoded in |10), 
gives a fidelity of 0.997. For larger numbers of qubits, the fi- 
delity scales favorably. For three and four qubits, scaling up 
the evolution time linearly with qubit number, we find fideli- 
ties of 0.989 and 0.990. 

The performance of the neutral-atom platform for AQC de- 
pends on a combination of practical and fundamental ques- 
tions. Here, photon scattering is the fundamental limitation on 
the fidelity. The size of the minimum gap for a certain level of 
dressing by the laser fields determines that time scale for im- 
plementing the algorithm and thus the probability of sponta- 
neous emission. The performance our protocol is constrained 
by the size of J due to Rydberg dressing, whose optimal value 
for a given laser power depends on the details of the atomic 
level structure. Here we found that for reasonable power and 
detuning, we could achieve J = 470 kHz and a fidelity of 
~ 0.99 in a proof-of-principle solution to a QUBO problem 
with 2-4 qubits. Modest increases in this coupling would al- 
low us to attain high-fidelty control with larger numbers of 
qubits. Alternatively, one could live with a slightly reduced 
fidelity (e.g., by increasing the detuning) at the expense of a 
modest increase in the number of times necessary to run the 
algorithm. This should allow us to explore the regime of 10-20 
of qubits, where interesting physics beyond classical simula- 
tion is accessible. The scalability this platform will depend 
on many factors, including a detailed study of many-body ef- 
fects beyond the pairwise interactions treated here lfT2l . Ul- 
timately, truly scaleable AQC will depend on our ability to 
harmonize the intrinsic robustness of the adiabatic approach 
with the the digitization of errors employed in the quantum 
circuit approach lfT5ll . 
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